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Abstract — We consider the problem of estimating the locations 
of a set of points in a /c-dimensional euclidean space given a subset 
of the pairwise distance measurements between the points. We 
focus on the case when some fraction of these measurements can 
be arbitrarily corrupted by large additive noise. Given that the 
problem is highly non-convex, we propose a simple semidefinite 
programming relaxation that can be efficiently solved using 
standard algorithms. We define a notion of non-contractibility 
and show that the relaxation gives the exact point locations when 
the underlying graph is non-contractible. The performance of the 
algorithm is evaluated on an experimental data set obtained from 
a network of 44 nodes in an indoor environment and is shown 
to be robust to non-line-of-sight errors. 

Keywords: Non-Line-of-Sight localization, semi-definite 
programming, robust matrix decomposition. 

I. Introduction 

The problem of localization has applications in many 
interesting areas such as cyber physical systems 
molecular biology nonlinear dimensionality reduction 
|3| etc. In all these applications, pairwise noisy distance 
measurements are obtained between subsets of points/nodes 
in some Euclidean space and it is required to estimate the 
locations of these points with high accuracies. Given the 
range of applications, the problem has received considerable 
attention in the last decade and is an active area of research. 
Even though the problem statement looks deceptively simple, 
the general case is shown to be NP-complete ^4J, i.e. finding 
a valid configuration of points satisfying a subset of pairwise 
distance measurements is not solvable in polynomial time. 
Significant research work is devoted to obtaining convex 
relaxations for the underlying optimization problem and 
obtaining conditions under which the problem can be solved 
in polynomial time. Semi-definite programming (SDP) 
relaxations have been proposed jSj, |[6j, jTj and the relaxation 
has been shown to obtain the exact solution for certain classes 
of graphs known as uniquely localizable graphs, examples of 
which include d-lateration and random geometric graphs. 

Most of the existing literature on convex relaxations focus 
on the case where the pairwise proximity measurements 
are either known noiselessly or known with some slight 
perturbations (i.e. line-of-sight (LOS) localization). However 
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(a) Sensor Localization. (b) Molecular conformation. 




(c) Manifold Learning. 



Figure 1. (a) Sensors placed in a field for monitoring the region, (b) 
Tryptophan, one of the 20 standard amino acids | 2 1. (c) (1) A high dimensional 
"swiss roll" data set that spans a lower dimensional hypersurface 1 8 1. (2) Data 
points within a local radius are connected. (3) The lower dimensional repre- 
sentation of the high dimensional surface preserving the relative geometry of 
the points. 

in practice, a significant fraction of the data could have 
outliers. For example, in the case of sensor localization, 
non-line-of-sight (NLOS) propagation introduces large errors 
in the distance measurements. Spin-diffusion phenomenon 
in NMR spectroscopy for protein molecular conformation is 
known to produce bias errors in the measurements |9|. Data 
in higher dimensions do not all lie on the manifold and there 
could be outliers that lie far away from the surface of the 
manifold. 

Traditional methods that do not account for the large 
bias errors can result in significant localization errors. The 
hardness of the problem arises from the fact that the bias 
errors are typically large and it is not known a priori which 



measurements are biased. There has been some work in the 
literature that addresses this case using convex relaxations 
i flQJ . However these relaxations assume that the locations of 
the errors are known which makes the problem considerably 
simpler to solve. Our focus here is to address the problem 
of localization in the presence of large unknown bias errors. 
For simplicity, we will keep the sensor localization problem in 
mind as we progress through the rest of the paper, though all 
the techniques can be easily adapted to other applications. We 
propose a semi-definite programming (SDP) relaxation for the 
case of NLOS localization and provide conditions under which 
the algorithm retrieves the node locations. The performance 
of the algorithm is evaluated through simulations and also 
validated on a data set collected through a real world indoor 
experimental setup of 44 nodes that is publicly available [llj . 
The proposed relaxation is shown to mitigate the effect of the 
NLOS errors to a large extent by obtaining accurate estimates 
of the node locations. 



II. Problem Setup 

Consider a static random placement of N nodes in a 
/c-dimensional euclidean space. The locations of these nodes 
are unknown and these nodes would hence be referred to as 
agents. We also have M special nodes whose locations are 
exactly known called as anchors. Nodes within a radius r of 
each other obtain relative distance measurements that could be 
corrupted by noise. Based on these distance measurements, 
the problem is that of determining the locations of the 
agents. Each of these measurements is modeled as either 
a LOS-dominated signal or a NLOS-dominated signal by 
choosing the observation noise to be drawn from a mixture 
of two distributions. 

Let Xi G denote the location of the ith agent. Let dij = 
\\xi — XjW denote the actual distance between the ith and jth 
node, dij denotes the corresponding distance measurement. We 
will model the distance measurement as df^ = df^ + bij + riij , 
where riij is a zero mean additive noise (e.g. gaussian) and 



bij is given by. 



bij = if LOS 
> otherwise. 



Under NLOS, bij is usually taken to be a positive random vari- 
able to model the statistics of multipath noise. The noise model 
is motivated by the time-of-arrival measurement modality for 
distance measurements. Since the primary multipath signal 
travels a longer distance than the true distance, the time delay 
of arrival for the NLOS path is larger than the LOS path and 
hence the NLOS bias can be modeled as additive. Similarly, 
the gaussian noise component is principally due to thermal 
noise in the receiver and so we take that term to be additive 
as well. We will assume that a fraction of the measurements 
are NLOS. We discuss a semi-definite programming approach 
to tackle this problem in the sections to follow. 



Algorithm 1 SDP relaxation for NLOS localization 

1: Input: distance measurements dij, G Sx U 
anchor locations {ak\^^i dimension k. 



minimize 

XX,B 



subject to 



Y X 



- d% + hj 



bij > 0, ij e Sx^Sa 



Y X 



y 0. 



3: Output: minimizer X. 



III. SDP Relaxation 

Given the problem statement, we can define a network 
{X, A^Ex^£a^ D), of N agents and M anchors, where X = 
[xf ; ^2^; x^] G R^^^ is the matrix of agent locations, 
A = [af ; ...;a^] G R^^^ is the matrix of anchor locations, 
£x is the set of edges between agents, £a is the set of edges 
between agents and anchors and D is the matrix of pairwise 
distances between nodes that belong to the edge set. An edge 
is present between two nodes in the network iff we have a 
distance measurement between them (e.g. if they are within a 
distance r of each other). We formulate the problem of esti- 
mating the node locations as an optimization problem. For this 
purpose, we will treat the biases bij's in the measurements, as 
parameters that need to be estimated. The maximum likelihood 
formulation for estimating the node locations as well as the 
biases gives us the following optimization problem. 



maximize 

X,B 

subject to 



bij > 0, iJ G £x ^£a' 



(1) 



The above optimization reduces to the following least 
squares minimization if we take n^j's to be i.i.d gaussian. 



minimize >^ 

X,B ^ 



^j\\2 



-dij + b,j 



( 



E 



V 



(2) 



subject to 



bij > 0, iJ e £x^£a 



where B = {bij} is the matrix of bias variables. The above 
minimization problem is highly non-convex and has multiple 
minima. We aim at obtaining a convex approximation to the 
above problem. 

The distance between the agents can be expressed as \ \xi — 



eJ^XX^ea where is a vector in 



pN 



with a 1 



and —1 at the zth and jth locations respectively and zeros 
elsewhere. Using the substitution Y = XX^, we get \\xi — 
^illl = ^Ij^^ij- The minimization problem ^ can thus be 
written as, 



Agent positions in 3D with the 
same graph structure and inter 
node distances 



minimize 

X,Y,B 



subject to 



[ejYci — 2eJ Xaj + aja^ 

hij > 0, ij e Ex^Sa 
Y = XX^, 



(3) 

where is a vector with a one at the ith position and zeros 
elsewhere. 

Since the constraint Y = XX^ is not convex, we relax 
it to F y XX^. Using Schur complements, we obtain 
the minimization problem as shown in Algorithm 1. The 
above semi-definite program can be efficiently solved using 
standard convex optimization techniques. One of the main 
advantages of this algorithm is that we do not need to input 
noise parameters to evaluate the optimization problem, which 
makes the algorithm useful in a practical setting. 

Let us develop an intuition for the relaxation. Without the 
relaxation, in problem ([3]), we were searching for a feasible 
set of points in the /c -dimensional space that minimize the 
cost function. The constraint that the points need to be in 
the /c-dimensional space essentially makes the problem non- 
convex and hard to solve. The condition Y >z XX^ relaxes 
the search for the points in any /c^ -dimensional space, such that 
k' < N . This can be seen as follows. If k' is the rank of Y , 
then Y h XX^ implies that there exists X' G IR^xC^'-^) 
^ X^ ^ 

can be thought of as the set of points in k^ dimensional 
space that are the minimizers of the optimization problem in 
Algorithm 1. Figure [2] illustrates this relaxation. The price that 
we pay for this simplification is that, once we get a solution 
in some k^ dimensional space, the projection onto the original 
/c-dimensional space could yield bad results. The rest of the 
paper explores conditions under which this relaxation actually 
gives a solution in the /c-dimensional space that would be close 
to the true node locations. 



A. Non-cooperative case 

Consider the simple case of a single agent and multiple 
anchors, i.e. N = 1 and M anchors. The minimization 
problem can be written as. 



such that Y = [XX'] 



X = [XX'] e 




Agent positions 
in 2D 



^ Anchor locations in 2D' 



Figure 2. Illustration of the SDP relaxation. We search for points in a higher 
dimensional space satisfying the same distance measurements and the graph 
structure as in the lower dimensional space. 



M 



minimize \^ (y — 2a[x -\- afai 

y,x,{bi} ^ V 

subject to bi > O^y i 



(4) 



y > x^x 

where x is the agent location to be estimated. 
The following theorem states our result. 

Theorem 1. For the case when the measurement noise is only 
due to bias errors i.e., rii = 0, the minimization problem 
gives an exact solution as long as the agent lies within the 
convex hull of two or more anchors with which it has LOS 
measurements. 

Proof: Let x be the actual node location, and suppose 
it is contained within the convex hull of two or more anchor 
nodes with which it has LOS measurements. By contradiction, 
let X be a solution to equation ^ such that x ^ x. Let W 
be the set of all i such that (if is a LOS measurement and 
X is contained in the convex hull of {a^,z G W}. Since the 
actual location of x is a solution, the optimal objective value 
is upper (and lower) bounded by 0. Thus it must be true that 
\\x - ai\\^ < (Jf = dD.Mi e W. For this condition to hold, 
{x — x^ai — x) > 0,Vz G W (consider 3 anchors with x in 
the middle. To be a distinct solution, x must simultaneously 
be closer to all 3 anchors than Xi's giving the inner product 
condition). But this is precisely the condition for the existence 
of a separating hyperplane between two convex sets, where the 
sets are the point x and the convex hull of the LOS anchors. 
If there exists a separating hyperplane between these two sets, 
then X is not in the convex hull of {a^,z G W}, which is a 
contradiction. So x must be the unique solution. ■ 

B. Cooperative case 

In this section we consider the cooperative case where 
we have internode distances between the agents. Recall the 
optimization formulation in Algorithm L Note that for every 
term in the cost function we have a positive bias term bij, 
that we wish to estimate and use it to offset any NLOS bias 
in the corresponding distance measurement. The optimization 
formulation does not assume any prior knowledge on the bias 
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Contracted 
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magnitude or whether a particular link is LOS or NLOS. 

Consider a given network of anchors and agents 
Ml = {X, A,Ex,£a^ D), with an underlying graph of 
connections between the nodes, where we have an edge 
between two nodes iff we are given a distance measurement 
between those two nodes. Suppose there exists another 
set of agent locations satisfying the same underlying 
graph topology as the original network A/i but with inter- 
node distances restricted to the edge set, shorter than the 
original network of nodes i.e., A/2 = {X^ A^Sv^Saj D), 
where Dij < Dij^ii^j) e £x ^ ^a- Figure pfa) is an 
example of such a network. In such a case, the optimization 
problem cannot distinguish whether the given set of distance 
measurements belonged to the original network, or whether 
the distance measurements were obtained from the second 
network corrupted by positively biased noise. Thus we can 
only hope to have a unique solution as long as the original 
network has a graph structure that rules out the existence 
of a different set of node locations satisfying the same 
graph structure and smaller inter-node distances. A fully 
connected graph is one such example. Figure [3jb) shows 
another example of such a network. Note that the node 
locations of the agents cannot be shifted without increasing 
at least one of the distances. Both these examples are in 
two-dimensions where we assumed that the nodes are also 
placed in two-dimensions (i.e. k = 2). However we saw that 
the SDP relaxation allows us to search in for node locations 
in a higher dimensional space. Hence we would need this 
property to hold in higher dimensions too. This motivates the 
definition of non-contractible networks. 




(a) Contractible network 



(b) Non-contractible network 



Figure 3. (a) Contractible networks. The dotted lines and corresponding 
nodes represent a contracted network with inter node distances smaller than 
that of the original network. (b) Non contractible network. See that the agent 
locations cannot be shifted without increasing at least one of the distances. 




3( (Proportion of NLOS readings) 



Figure 4. Error as a function of the fraction of NLOS measurements under 
the different formulations. The NLOS bias was uniform on [0,6], the noise 
was zero-mean Gaussian with ctlos — 0.02, and error values are averaged 
over 10 trials each. 



Definition: A network (X, A, Ex^^a-, D) is said to be non- 
contractible if there exists a unique X G R^^^ satisfying 
the location constraints imposed by D and there exists no 
x'j e M^, j = 1, 2, .., N,\/h>k such that 

||(a,;0)-4lP < 4 V {iJ)eSA 
< 4 V {iJ)eSx 
Xj ^ {xj;0) for some j G {1, 2, 3, TV} 

Figure [3] shows examples of networks that are contractible 
and non-contractible in R^. Note that even if one node 
lies outside the convex hull of the anchors, the network is 
contractible. However the converse does not hold. 

The following theorem establishes that non-contractibility 
is a necessary and sufficient condition for the SDP relaxation 
to provide an exact and unique solution. 

Theorem 2. For the case when the measurement noise is only 
due to bias errors i.e. Uij = 0, 

• Algorithm 1 gives the exact solution if the underlying LOS 
network of nodes is non-contractible. 

• If the max rank solution of the optimization problem in 
Algorithm 1 is k, then the solution is exact. 



Proof: The intuition for the proof relies on the arguments 
made for the definition of non-contractibility. The sufficiency 
of the non-contractibility condition is proved using arguments 
of contradiction. See the appendix for a formal proof which 
is on the lines of ||5|. 

■ 

Remark: The notion of non-contractibility is a general- 
ization of the uniquely localizable condition |5 | for NLOS 
localization. However, the question of what networks are non- 
contractible in practice is still an open question. There has 
been recent progress in the LOS case in characterizing the 
classes of networks that are uniquely localizable using graph 
rigidity theory |^. Analogously, one could hope to utilize 
results from tensegrity theory |12| , to characterize networks 
that are globally stable and non-contractible. This is a future 
direction of research. 

IV. Validation 

A. Simulation 

Simulations were carried out for the proposed SDP 
relaxation and compared against the scheme proposed by 
Biswas et. al. |5 |. We chose this baseline system to compare in 
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Figure 5. Per-node error as a function of communication radius on the data 
set obtained by Patwari et al |,13j . 
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Figure 7. Node location estimates for the experimental data set using the 
scheme by Biswas et al. 0. 
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Figure 6. Node 
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location estimates for the experimental data set using our 



order to highlight the gains obtained by expHcitly accounting 
for NLOS in the measurements. The SDP relaxation in |[5| 
is robust to small percentages of NLOS noise but fails for 
larger percentages given that they do not explicitly account 
for NLOS. The setup consisted of 50 nodes randomly placed 
in a [—1,1] X [—1,1] grid with 15 anchors. Anchors were 
placed uniformly along the upper and lower boundaries of 
the grid. The radius of connectivity was taken as r = 1.5, 
and the gaussian noise standard deviation was set to 0.02. 
The NLOS noise was taken to be uniform over the interval 
[0,6]. Results were obtained after averaging over 10 trials. 

Figure]?] shows the variation of the error in the node location 
as a function of the fraction of NLOS measurements for the 
relaxation. The error metric is the per node average mean 
squared error in the node location estimate i.e. \\x — x\\/\^. 
It can be seen that the proposed SDP relaxation is robust to 
a large fraction of NLOS errors. 

B. Experimental Results 

The SDP relaxation was tested against data obtained from 
real world experiments. We use the data set published by 
Patwari et. al. |13|. The experiments were conducted in an 



indoor office environment by simulating 44 node locations 
using a transmitter and a receiver and obtaining pairwise 
time-of-arrival measurements. The environment had a lot 
of scatterers and nearly all the distance estimates have 
strong NLOS biases in them. The authors in [ |13| , artificially 
subtracted out the NLOS biases and obtained the node 
location estimates using LOS algorithms. They report a 
per node localization error of 1.26m after utilizing all the 
measurements i.e. a fully connected network. We use the 
same data set without subtracting out the NLOS biases and 
show that the SDP relaxation performs very well even in this 
setup. The node location estimates are obtained using only 
a subset of measurements, based on an arbitrarily defined 
radius of communication. The SDP relaxation does not need 
any inputs regarding the noise parameters which is required 
by most of the existing algorithms in the literature. 

Figure [5] shows the per node error as a function of the 
communication radius. The per node error is seen to be less 
than a meter even when the NLOS biases are included and 
only a subset of the measurements are used. However the 
performance declines as the radius becomes too large or too 
small. The reasoning is that for small radii, each node may 
not have enough information to fully resolve its position and 
the network could become non-contractible. For large radii, 
the multipath bias is more likely to be large, causing errors 
which can be reasoned as follows. The measurement model 
that we considered has the bias added to the square of the 
distances. However in practice the bias would be added to the 
actual distance value. Hence when we square the distances 
and apply our measurement model, the bias values would be 
a function of the distance and hence would increase as the 
distances increase. This partially explains the reason why the 
performance degrades slightly as we increase the radius of 
connectivity. Figure [6] and Figure |7] show the estimated and 
true node locations for the experimental data set obtained from 
the SDP relaxations. 

V. Conclusion 

In this paper, we considered the problem of estimating the 
locations of a set of points in a /c-dimensional euclidean space 



given pairwise distance measurements amongst the points. We 
considered the case when some fraction of the measurements 
could be arbitrarily corrupted, proposing a SDP relaxation 
formulation for the problem. The relaxation was shown to 
give the exact solution when the underlying graph is non- 
contractible. Simulations and testing on real-world data shows 
that the algorithm is quite robust to NLOS noise. Future 
work involves strengthening the theoretical results in order 
to provide meaningful guarantees in practice. It would be 
interesting to characterize the classes of graphs that are non- 
contractible. We also need to theoretically prove that the SDP 
relaxation is robust to thermal noise. 
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We will now show that if the underlying LOS network 
is non-contractible then the solution is unique. Assume that 
each node has at least one LOS measurement. We will prove 
uniqueness by contradiction. Suppose that the network is 
non-contractible and the solution is not exact. Then by the 
definition of non-contractibility and by the non-uniqueness of 
solution X, there exists a solution X' of rank higher than k 
that satisfies Y y X'^ X' . X' satisfies 



Y - X^X = X'^X'. 



Define 



We then have 



X 



X 
X' 



X^X 



■X'^X' 



x^x 

= Y. 

Let Slos be the edge set of the LOS links. 

0, 



■2 

elYe, 



4 



< 
< 



4 



Vi,i G £los^£x, 
eSLOsnSx- 



We also have, 



ejYe, 
2[aj0] 



Y X 
h 

2aJXe. 
X 
X' 



T 

■ Ci£ (l£ 



Ikll' 



(a^;0)-x,|p 
for 



< 
< 
< 
< 



4 
4 

£los n £a 



'2 



Thus, we just showed that the underlying network is con- 
tractible which is a contradiction to our initial assumption that 
the network is non-contractible. 



Appendix A 
Proof of Theorem[2] 

Lemma 1. If the max rank solution of the optimization 
problem in Algorithm 1 is k, then it is unique. 



Proof: We will prove this by contradiction. Suppose 
(Xi,Yi,5i) and (X2, 12,^2) are two solutions, then for 
any /3 < 1 one can easily check that /3(Xi,Yi,5i) + (1 — 
/3)(X2, 12, ^2) is also a solution. Since the max rank is k, we 
should have that. 



f3Y^ + (1 - /3)Y2 
pX, + (1 - p)X2 

be rank k. This gives us 



Ik 



(1 



■mi 



(/3Xi + (1 - /3)X2f{/3Xi + (1 - /3)X2). 

T ■ 



PYi + (1 - P)Y2 

where Y\ = X^Xi and Y2 = Xj X2- Rearranging terms, we 

(2 



get ||Xi — = which implies Xi = X2. 



